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^ ! Abstract 

^ ^ ' A new algorithm for calculating the spin tune and the n-axis for circular accelerators is 

^' presented. The method resembles the one employed in the existing program code SODOM 

^ • in that one-turn numerical spin maps at equally spaced orbit angle variables are used but 



it is more efficient than the latter. Furthermore, it is applicable at large openning angles 
of the n-axis, whereas the existing SODOM only converges for small angles. 



^ On leave of absence from the National Accelerator Research Organization, Japan (KEK). 



1 Introduction 

The spin tune, the spin precession frequency divided by the orbit revolution frequency, is 
an important parameter in the description of spin motion in circular accelerators. When 
a particle is on the closed orbit, the definition of the spin tune is obvious; it is the 
spin precession angle over one turn divided by 27r. However, when orbit oscillations are 
involved, the definition of the spin tune becomes more complicated. One needs the concept 
of the n-axis which was first introduced by Derbenev and Kondratenko[jl| for radiative 
polarization phenomena in electron storage rings. 

We assume that we have complete knowledge about the orbit motion, i.e. that we 
know the action and angle variables, J = {Jx, Jy, Jz) and = {(px, <Pyj (pz), corresponding 
to the three degrees of freedom of the orbit motion, which can be in general nonlinear. 

A particle with initial coordinates (J, (j)) at a machine azimuth 6 executes orbit oscil- 
lations and comes to (J, (/> + fi) after one turn {9^9 + 27i), where fj, = {fix, fiy, fiz) is the 
orbit tune, u, times 2n. The spin motion over one turn can in general be expressed by a 
3x3 rotation matrix R{J, 0, 9). Obviously, it is a periodic function of 9 and (/) with period 
27r. On the next turn the rotation is expressed by R{J, (f) + fj,,9 + 2tt) = R{J, (f> + fj,,9) 
which differs from R{J, (f), 9) unless the orbit tunes are integers. 

A particle on the closed orbit sees the same rotation Ro{9) for every turn. Rq{9) 
has eigenvalues 1 and e^*^^° and the spin tune Vgo is /iso/27r. One can show that /Xso is 
independent of 9. The eigenvector belonging to the eigenvalue 1 is denoted by no, i.e., 
RqUq = tiq. It depends only on 9. A spin parallel to no{9) remains unchanged after one 
turn, and all other spins attached to closed orbit trajectories precess by the angle figo 
around no during one turn. 

The vector n is a generalization of no for particles off the closed orbit. It is a function 
of (J, <f), 9) periodic in (f> and 9 and satisfies 

R{J,<l>,9)n{J,(l),9)=n{J,(l> + ti,9). (1) 

When J = 0, n reduces to no. To define the spin tune for nonzero J, we need two more 
vectors Wi and U2 which form an orthonormal basis together with n. They are functions 
of (J, (f), 9) and periodic in <f> and 9 like n. The spin tune is defined as the precession angle 
in the frame {ui,U2,n) divided by 27r. 

The concept of the vector n has been playing an important role in the description and 
calculation of radiative polarization in electron/positron storage rings since|]l]]. Recently, 
it has also turned out to be useful for proton rings p|. 

To calculate the vector n several algorithms have been invented. S. Mane0 developed 
a computer code SMILE using a perturbation expansion with respect to the orbit action 
variable. The present author suggested a perturbation algorithm using Lie algebra]^ and 
Eidelmann and YakimenkoQ coded a program SPINLIE with (low order) orbit nonlin- 
earity. Balandin, Golubeva and Barber]^ also wrote a Lie Algebra code. 

The present author considered another method which does not employ a perturbation 
expansion and wrote a program SODOM0]. Heinemann and HoffstaetterP] use tracking 



and 'stroboscopic averaging' in the code SPRINT. The programs SODOM, SPRINT and 
additionally compute the spin tune. 

The new method which we are going to describe is based on SODOM. 

We shall briefly summarize the SODOM algorithm in the next section and describe 
the new method in Sec. 3. 

2 The SODOM Algorithm 

Let us first briefly summarize the algorithm employed in SODOM. (See Sec. 3 of [0.) 
Denote the one-turn SU2 spin transport map starting at a fixed prechosen azimuth 6q for 
particles with initial orbital phase (j) by M(0) and the spinor representing the n-axis at 9o 
by '?/'(0). (Here, we simply write ip^cj)) instead of i'+{(j>) ^- We also omit the arguments 
J and 6q since we shall deal with one set of J and consider the one-turn map from the 
origin ^o only.) The fact that n is 'invariant' means 

M(0)^(0)=e"^^(^)/27A(0 + /^), (2) 

where f (0) is a real periodic function. Once a solution {iIj{^),v {</))) is obtained, we solve 
the equation 



and define 



Then, \E'(0) satisfies 



v{(f>) + u{(f) + fj.) - u{(f>) = fis (3) 



^(0) = e'"^'^)/^!^) (4) 
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M{(l>)^{(j>) = e-'^=/2^(0 + /i), (5) 

where /i^ is the spin tune times 2tt. The Ui 2 axes are represented by a spinor 

g-iv/2^ + e^^/^§*l (6) 

where we define the operation ^ as 

§ = 2^2^*, (7) 

which was denoted by \E'_ in 0. Note that \E' = — \E' and \E'^\E' = 0. The three spinors, "^o, 
'^n/2, ^; represent the three vectors Ui, U2, n. The phase of \E' is irrelevant for defining n 
but it is important for Ui and U2- 

The original SODOM algorithm parametrizes ip as 
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The SU2 matrix M(0) can be parametrized by two complex functions /(0) and g{^) as 
Then, one gets an equation for ^: 

^* wcw + 9{<t>)a<i> +ti) = f{<i>) - r{<t>)a<t>)a<i> + /^)- (lo) 

By expanding /(0), (7(0), and (^(0) into Fourier series hke I] /me*"*'*^, we get a nonhnear 
equation for Cm- 

A key component of SODOM is the calculation of the Fourier coefficients fm and 
Qm from the tracking data over one turn for several particles having the same J but 
equally-spaced (0 < < 27r). 

The parametrization (||) is good only when (,{4') is small. Because of its up-down 
asymmetric form, many more Fourier terms are needed than required by the physics, 
when C(0) is large.0 Also, the iterative method of solving the nonlinear equation easily 
fails when C, is large. 

3 The Matrix Eigenvalue Method 

The new algorithm is much simpler and involves solving eq.(^) directly rather than eq. (p!OD . 
By expanding M{<f)) (actually the functions /(0) and g{(l))) and \l/(0) into Fourier series 
as 

M(0) = Y. M^e^-^ vl/(0) = ^ vl/^e^-"^ (11) 

171 TTl 

eq.(^) can be written as 



e— '^ Y: Mm-m'^m' = e-^'-^/'^m. (12) 



This is simply a matrix eigenvalue equation. Thus, the spin tune comes out as an eigen- 
value. 

However, obviously, eq.(0) has many eigenvalues. Which one gives the spin tune? 
What do the other eiganvalues and eigenvectors mean? In order to answer these questions 
let us return to eq.(|D and examine it as an eigenvalue system 

M(0)^(0) = A^(0 + //) (13) 

Note that this is not a simple 2x2 algebraic equation because of the (f) + /J,. 

Before going further we have to think about subtle problems associated with the '2-to- 
r correspondence between SU2 and 803. Note that we use 2-component spinors and SU2 



^ For example, ^ — (cos 0, sin 0) is a mild function but leads to C, = tan0 which is hard to Fourier- 
expand. 



matrices instead of 3-vectors and 803 matrices to achieve computational speed and to 
minimize storage but not because the particles have spin %/2. The classical spin motion 
can be completely described by 3-vectors and S03 matrices. 

What does the periodicity of a spinor with respect to mean? The physical object is 
the 3- vector ^'^a^ = n rather than the spinor \1/. In this sense a complex phase factor in 
\1/ is irrelevant. However, a complex phase factor is still relevant when one constructs the 
Ui and U2 axes from \1/ via \1/^. 

On the other hand, a sign change of \I' does not cause a change of n = \1'V\1' nor a 
change of Ui and U2 defined by \l/^<7\l/(p. 

Thus, as the periodicity condition for \1' with respect to (pj (one of the orbit angle 
variables), we have to allow both '^{(j)j + 27r) = \I^(0j) and '^{(j)j + 2-n) = -^(0^). Then 
with 3 degrees of freedom for orbit motion, we have 8 types of solutions \l/(0) differing by 
their sign change behaviour under the transformation (pj -^ (pj + 2tt. In Fourier expansion 
language, this means that \l/ can be expanded as 

^(0) = e^'""-^/^ ^ ^^e'"'t> (14) 

m 

where m° is a set of three integers each of which is either or 1 . 

We now define the scalar product of two arbitrary spinors \E'i and \l/2 by 

1 f^'^ J- J- 

(Vl>i,vl>2)^— - / Vl/t(0)vl/2(0)# = 5w>o^<„,*2,r. (15) 

(4vr)^ Jo ' ' m 

Obviously, solutions of different types in eq. (|l^ are always orthogonal. In the following 
we consider the solutions of eq. ([T3|) which are 'periodic' and smooth in 0. 



Now, let us list a few lemmas. 

[a] |A| = 1 

[b] (^i,^2) = 0z/Ai^A2 

From the unitarity of M(0) we get 

vI/J(0)M/,(0) = vI/J(0)M(0)tM(0)M/,(0) = [Af(0)vI/,(0)]Uf(0)vI/,(0) 

Integrating over (j) and using the definition (p!5|), we get [a] for i = j and [b] for Aj 7^ Xj. 
Note that [b] does not imply ^1(^)^^2(0) = for Aj ^ Xj. 

[c] 1^(0)1 is independent of (j) (and can be normalized to unity). 

The unitarity condition |\I/(0)| = |\l/(0 + ^)|, together with the smoothness of \l/(0), and 
the non-commensurability of fi are enough to guarantee [c]. 

[d] // (A, ^(^)) is a solution, so is (A*, ^(0)) 

Take the complex conjugate of eq.(0) and use a2M*a2 = M. If \1/ corresponds to 
n, then '^ corresponds to ~n and the spin tune changes sign. (Since a2(Ta2 = —c*, 
$t£r§ = -vpto-vl/.) Note that not only (§, \&) = but also §"f^ = at every (f>. 



[e] If \ is an eigenvalue, then so is Ae*'"''*/^, where m is a set of any integers. 
Multiply eq.(|T3p by e-*"*-^/^ ^nd define §(0) = e''"''^^^'^ {(f>) . Then 

M{<l))^{(l)) = Ae-''"-^/2^(0 + /i) = Ae''"-''/2^(0 + //) 

Thus, ^ is an eigenvector belonging to the eigenvalue Ae*"*'*/^. 

This gives an ambiguity in the spin tune: fig -^ fis + T^fJ,. However, all the eigenvalues 
of the form Ae*"*''/^ give the same vector n = \E' V\l' = \1' V\E'. The Ui^2 axes corresponding 
to ^ are 
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g-*V'/2g-im-<^/2^ _|_ ^iip/2^im<t>/2--qj* 



^ ip+m-^ 



Thus, the new Ui 2 axes rotate by m-0 with respect to the original ones. 

From the lemmas above, we know that once a solution (A, \1') is found, we can construct 
infinitely many solutions of the form (Ae*"*''/^, e-*'"'^/^^) and (A*e-*'"''/2^ e*'"'^/^^), and 
that they all correspond to the same vector n or —n. 

A natural question is then 'are there any other eigenvalues?'. The answer is 'No': 

[f] // A is an eigenvalue, all other eigenvalues are either Ae*"*'*/^ or _x*e^*'"''/^ 

If (Ai, ^1) and (A2, ^2) are solutions, a(<^) = ^5(<^)^i(<^) = (M*2)"^M^i = A^Aia(<^+/x). 
From the periodicity and smoothness of a{<j)) and the non-commensurability of// one finds 
either that [1] a(0) = e*"*^ and AgAi = e~*°''^, a being a constant 3-vector, or that [2] 
a = (\E'i and \E'2 are locally orthogonal). In the case [1] a must be of the form m/2 from 
the periodicity requirement, where m is a set of three integers. Therefore, A2 = Aie*"***/^. 
In the case [2], examine ^^2 in place of \E'2. Then we get either A2 = AJe"*'"'^/^ or 
§tv|>^ = 0. However, if both \E'^\E'i and ^l^^i vanish, then \[^i = because \E'2 and §2 
are orthogonal. Therefore ^^2^1 = cannot be the case. Thus, the cases [1] and [2] 
correspond to A2 = Aie*"*''*/^ and A];e^*'"''/^, respectively. 

Let us consider the spin tune Ug = /is/27r. It is obtained from the definition A = 
g-«Ms/2 _ g-JTi/s^ Prom the above arguments we find that if the set [us,n,Ui + iu2] is a 
solution, then [ug — to-i/, n, e~*"*''^(ui + iu2)] and [— z/^ — mv, — n, —e^'^'^iui + iU2)] are 
also solutions. Thus, the spin tune has ambiguities up to a multiple of the orbit tunes 
and up to a sign. The latter is related to the choice of sign of n. 

When obtaining z/^ from A, one finds an ambiguity only up to an even integer rather 
than up to an integer. At first sight this is puzzling but it is also due to the '2-to-l' 
correspondence between SU2 and 803. Obviously, 

[g] //\E' is an eigenvector of M with eigenvalue \, it is also an eigenvector of —M 
with eigenvalue —A. 

Since M and —M represent the same SOS rotation, we have also to include the solutions 
to —M. However, —M has exactly the same eigenvectors as +M (therefore the same 
Ui,U2,n) with spin tunes z/^ shifted by one. This solves the above puzzle. Thus, we can 
define the spin tune in the interval [0,1) or (-0.5,0.5] and, if the sign of n is irrelevant, we 
can reduce the interval to [0,0.5]. 
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Thus, we have found that we only need one of the sets of eigenvector and eigenvalue of 
eq.(^). All others can be constructed from this. Therefore one can Fourier expand "$ as in 
sq-dUD rather than as in the general from (0). (Note, however, that one will find tunes 
of the form ±Us + 2mv although odd-multiple solutions can be easily reconstructed.) 

Let us briefly discuss the degeneracy. Within the eigenvalue group Ae*"*'*/^ of the 
same sign of n, a degeneracy is possible when mv is an even integer, i.e., when the orbit 
motion is in resonance, which we are not interested in. We may assume this is not the 
case. 

On the other hand, a degeneracy between solutions of different signs of n corresponds 
to a spin-orbit resonance. When the two solutions (A, \E') and (A*, \E') degenerate (A = A*), 
the spin tune becomes an integer. Taking into account the ambiguity of spin tune, this is 
equivalent to the relation z/^ = mu + integer. 

4 Choice of the Spin Tune 

We have shown in the previous section that there are many eigenvalues (spin tunes) 
representing the same vector n and different {ui,U2) axes. Now, we must finally decide 
which eigenvalue to choose for the spin tune. Theoretically speaking there is no reason 
to choose one particular value. As pointed out in spin tune is intrinsically ambiguous 
up to a multiple of the orbit tunes. The choice of the spin tune, which is equivalent to a 
choice of {ui,U2), is a matter of convention. 

In practice, however, a solution is not desirable if (111,1*2) is a strong function of 0. 
When one solves the equation by Fourier expansion, the most natural choice is to take 
the solution having the largest zero-Fourier harmonic \E'm=(o,o,o) • 

If one plots all the eigenvalues as a function of any parameter (beam energy, betatron 
amplitude, etc), one will find continuous curves. If one plots the spin tune selected as 
just described as a function of these parameters one may occasionally find a jump of spin 
tune although the whole spectrum content is continuous. 

Let us give an example from a test calculation. The test ring consists of 100 FODO 
cells, each of which has two thin-lens quadrupole magnets and two bending magnets 
filling the entire space between quadrupoles. The focusing effect of the bending magnets 
is ignored. In order to avoid a too high symmetry of the orbit motion, an artificial phase 
advance of 90 degrees in both horizontal and vertical planes is introduced at one point 
in the ring. The tunes are u^ = 15.3827 and Uy = 25.6482. Only the vertical betatron 
oscillation is excited. The beam energy is so chosen that z/^o = 7a = 1520.72. 

Fig.l shows the eigenvalue (spin tune) spectrum as a function of the betatron action 
Jy. Only those with small m are plotted. The points linked by a solid line correspond to 
the spin tune selected by the criterion mentioned above. As one can see, each eigenvalue 
is a continuous function of Jy (A few curves appear broken because not all the eigenvalues 
are plotted.) but the selected tune shows a jump at Jy ~ 0.7 x 10~^m-rad. The spin tune 
before and after the jump, Ugi and 1/^2 satisfies z/^i + 1^32 = —'^^y + integer. The dashed 



line (with the same scale) is the upper limit of polarization, i.e., 
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The minimum of Pum coincides with the point of the spin tune jump. 
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Figure 1: An example of a spin tune spectrum as a function of betatron action. The 
dashed line is the upper limit of polarization Pum- 

We have compared the results of our program with SPRINT for the amplitude depen- 
dence of the spin tune in the HERA ring. The agreement of the n axis and the spin tune 
was excellent. Not only the occurrence of spin tune jumps but also their location agree, 
which means that taking the largest zero-harmonic and the stroboscopic averaging are 
almost equivalent. 

5 Truncation of Fourier Series 



In numerical calculations one has to truncate the Fourier expansion. There are a few 
problems associated with the truncation. 



When A^ values of are used (we deal with one degree of freedom for illustration. 
The extension to 3 degrees of freedom is obvious.), the range of the harmonics should 
be — M < m < M {N = 2m + 1).0 For a discrete Fourier transform the range can also 
beO<'m<A^ — 1 (as in standard FFT routines), but this choice is not good when 
other values of (p are needed (for example when calculating n for arbitary values after the 
problem is solved). A^ must be large enough to ensure that the Fourier components Mm 
(actually fm and gm) are small enough outside the region [— m, m].0 

The matrix e~*™'^Mm-m' in eq.(0) is then a (2m+1) x (2a/+1) matrix (each element is 
a 2 X 2 matrix and we are dealing with one degree of freedom.). One finds that the diagonal 
elements (m = m') are normally large and that the elements with large \m — m'\ are small. 
The elements in the upper-right and lower- left triangles {\m — m'\ > m) are exactly zero 
because they require the harmonics outside [— m, m]. Owing to this truncation, the matrix 
does not exactly satisfy the lemmas in the previous section even if A^ is very large. (For 
example, in the first row {m = —m) even the first harmonic m — m' = lis lost because m' 
would be —M — 1.) Although the solution with the largest zero-harmonic is not affected 
much by this truncation, it is not easy to confirm the accuracy of the solutions. 

On the other hand, one can fill the upper-right and lower-left triangles by treating the 
harmonics in a cyclic manner as in a discrete Fourier transformation (i.e. one identifies 
the {m + l)-th harmonic with the (— M)-th.). With this prescription the truncated matrix 
becomes exactly unitary even if N is not large enough. The solution with the largest 
zero-harmonic does not change much. The appearance of eigenvalues with modulus far 
from unity then means that the eigenvalue solver is not accurate. 

When one adopts the cyclic use of the harmonics, the lemmas [a], [b], [d] and [g] hold 
exactly apart from round off errors, but [c] and [e] (and accordingly [f]) become inaccurate. 

6 Conclusion 

We have shown that the spin tune can be obtained as an eigenvalue of a matrix which is 
created from the one-turn maps calculated by particle tracking. The method is applicable 
to any system with linear or nonlinear orbit motion as long as the orbit action variables 
exist. The convergence is much better than with perturbation methods and the previous 
SODOM algorithm. The computation is very fast because it makes full use of the fact 
that the spin motion is linear and that we know the orbit tunes. 
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